function ans=Energy1(T)
    global Tstar;
    global ku;
    ans=ku*abs(Tstar-T);
end